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Multi-component generalizations of the CH equation: 
Geometrical Aspects, Peakons and Numerical Examples 
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Abstract 

The Lax pair formulation of the two-component Camassa-Holm equation (CH2) is generalized to produce an 
integrable multi-component family, CH(n,k), of equations with n components and 1 < < n velocities. All 
of the members of the CH(n,k) family show fluid-dynamics properties with coherent solitons following particle 
characteristics. We determine their Lie-Poisson Hamiltonian structures and give numerical examples of their 
soliton solution behaviour. We concentrate on the CH(2,k) family with one or two velocities, including the 
CH(2,-1) equation in the Dym position of the CH2 hierarchy. A brief discussion of the CH(3,1) system reveals 
the underlying graded Lie-algebraic structure of the Hamiltonian formulation for CH(n,k) when n > 3. 

1 Introduction 

Purpose. 

This paper develops the Lax formulation for an integrable family of equations that contains the multi-component 
generalizations CH(n,k) of the CH equation with n momentum components convected by 1 < \k\ < n veloci- 
ties. In the CH(n,k) family, CH is designated by CH(1,1) and CH2 by CH(2,1). We also consider the Lie-Poisson 
Hamiltonian properties of several other members of the CH(n,k) family, particularly CH(2,1), CH(2,2), CH(2,-1), 
CH(3,1) and CH(3,2). The CH(2,-1) system may be regarded as a two-component generalization for CH of the Dym 
equation in the CH hierarchy. All the CH(n,k) equations for a given value of n share the same spectral problem 
and therefore belong to the same integrable hierarchy. As might be expected, higher-order powers of the spectral 
parameter appear in the Lax formulation of CH(n,k) for greater values of n. When compared to rigid rotations, the 
CH(2,1) shallow water system recovers the heavy top equations and the CH(2,2) system recovers the equations for 
a rigid body in a potential field. The CH(3,1) system reveals the underlying graded Lie-algebraic structure of the 
Hamiltonian formulation for CH(n,k). Examples of numerical solutions of these integrable systems of equations 
illustrate their interesting dynamical properties, in which soliton trains emerge from spatially confined initial con- 
ditions and interact with each other in a variety of different ways for the various CH(n,k) systems investigated here. 
Many open problems arise and we attempt to sketch some of the opportunities for future research in the conclusion 
section. 

1.1 Brief review of the Camassa-Holm (CH) equation 

The CH equation EJ|71 

ut - u xxt + 2uju x + 3uu x - 2u x u xx - uu xxx = 0, (1.1) 

governs the evolution of the function u(x, f) : R x 1 I, interpreted as a shallow water fluid velocity. When the 
linear dispersion parameter to G M vanishes (lo = 0) the CH equation ( |1.1| > admits peakon solutions. Peakons are 
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nonanalytic solitons that superpose as 

N 

u{x,t) = iX>«(*) e " l<B " ? " (t)l ' (1.2) 



a=i 



for sets {p{t)} and {q(t)} satisfying a system of completely integrable canonical Hamiltonian equations. 
For the CH equation ( |1.1[ ) describes the unidirectional propagation of dispersive shallow water waves over 

a flat bottom at one order higher than KdV in the standard asymptotic expansion of the Euler fluid equations with a 
free surface in a certain Galilean frame (mean wave speed) ||6j|7j[T6l[T7j|39l|40l. It also describes axially symmetric 
waves in a hyperelastic rod lTT5l . The inverse scattering problem for CH is treated, e.g., in (4l[TT][T2l[T3l. For a 
brief history of the developments and results about the CH equation of relevance in the present paper, one may 
consult the recent review 11261 and references therein. 

The scope of the variety of mathematical interpretations of CH may be gleaned by rewriting it in various equivalent 
forms, many of which have been discovered several times before and each of which can be a point of departure for 
further investigation. For example, CH may be treated variously as: a fluid motion equation; a mathematical model 
of shallow water wave breaking; a vanishing Lie derivative, describing invariance of a 1-form density under the flow 
of a vector field related to it by inversion of the Helmholtz operator; a nonlocal characteristic equation; an Euler- 
Poincare equation describing geodesic motion on the diffeomorphism group with respect to the metric defined by 
the H 1 norm on the tangent space of vector fields; a Lie-Poisson Hamiltonian system describing coadjoint motion 
on the Bott-Virasoro Lie group; a bi-Hamiltonian system; a compatibility equation for a linear system of two 
equations in a Lax pair, etc. 

Our own point of departure and primary emphasis in this paper is in treating CH as a member of a family of inte- 
grable evolutionary equations associated with a certain class of energy-dependent isospectral eigenvalue problems 
of polynomial order. We shall also discuss the geometrical interpretations and varieties of solution behaviour of 
several members of the hierarchy. To establish notation, we begin by writing CH in several of its various forms and 
discussing its properties that are relevant here. Most of these properties were already established in ||6ll2l- 

CH as a fluid motion equation. CH may be written in the form of a fluid motion equation, as 

u t + uu x = -P x , (1.3) 

with pressure P given by the convolution 

P = K * (u 2 + \u 2 x + uju) with kernel K (x, y) = | exp(— \x — y\). (1.4) 
The kernel K is the Green's function for the ID Helmholtz operator, (1 — d 2 ), and is also the shape of the peakon 



profile in equation ( 1.2 1. 



Lie derivative, or characteristic form of CH. Upon introducing a momentum variable 

m = u - u xx = (1 - dl)u, (1.5) 

CH expresses the vanishing Lie derivative condition for invariance of its momentum, as a 1-form density, along 
characteristics of an associated velocity vector field, 

(dt + C u ) {{m + uj) dx 2 ) = (^rnt + um x + 2(m + oo)u x ^j dx 2 = , with u(x,t) = K*m. (1.6) 

This form of CH may be interpreted as the the condition for the 1-form density (m + uj) dx 2 to be preserved (frozen 
in) under the flow of a characteristic velocity dx/dt = u(x(t),t), namely, 

— l(m + uj)dx 2 j=0, along — = u(x(t),t) = K * m. (1.7) 

(JjL \ / CvL 
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In the parlance of fluids, equation ( 1.6 1 is the Eulerian form of the invariance, while equation ( |1.7[ ) is its equivalent 



Lagrangian form. As with the evolution of vorticity according to Euler's equations for incompressible fluid motion, 
the relation between the velocity of the flow and the property it carries is nonlocal. For the Euler fluid equations, 
the velocity carries the vorticity, related to the velocity by the convolution in the Biot-Savart law, which inverts the 



curl operation. For CH, the momentum is related to the velocity by convolution with the kernel K in ( L4l, which 



inverts the Helmholtz operator. When expressed in terms of the momentum, the peakon velocity solution ( 1.2) of 



dispersionless CH for w = becomes a sum over delta functions, supported on a set of points moving on the real 



line. That is, the peakon velocity solution ( 1.2 1 implies 

N 



m(x,t) = Y J Pa(t)S{x-q a (t)), (1.8) 



a=l 



because of the relation (1 — d 2 ) (ie ' x y ') = 5(x — v) for the kernel K in ( 1.4 1. As shown in 11271 . the CH peakon 



solution ( 1.8 1 is geometrically the cotangent-lift momentum map for the left action of the diffeomorphisms Diff (R) 



on a set of N points on the real line. 

The Euler-Poincare and Lie-Poisson properties of CH The geometric properties of CH in \IA\ follow from 
Hamilton's principle with a Lagrangian l(u) : X(M) — > E by using the Euler-Poincare theory ll28l . Its associated 
Lie-Poisson Hamiltonian formulation in terms of m £ I* (R) then emerges from a Legendre transformation. Here, 
the velocity vector field u G X(R) ~ TDiff (R)/Diff (R) in the Lagrangian l(u) is right-invariant under Diff(R), 
the diffeomorphisms of the real line. The Lagrangian for CH is l(u) = |||it|| m, the H 1 norm on these vector fields; 
which provides the metric for the interpretation of CH solutions as geodesic motion on Diff (R). The variational 
derivative 5l/5u G 3E*(R), in the space of real-valued 1-form densities dual to X{M) yields the CH momentum, 

51 

m = — = u-u xx . 
ou 

The Lie-Poisson Hamiltonian formulation of CH follows by Legendre transforming the Euler-Poincare equation 
for right-invariant vector fields, as 

d 51 1Jk 51 _ , 5h 

— — = -ad n — m t = -(d x m + md x )— (1.9) 

at 5u 5u dm 



with 



him) = (m, u) — Liu) , - — = u, — = V = m — — , 

dm ou 5u 



where ( • , • ) : X* (K) x X(M) — > M denotes L 2 pairing on the real line. For more detail of the Lie-Poisson 
structure of CH, see 11281 . For its interpretation as coadjoint motion on the Bott-Virasoro Lie group, see [49j. The 



Euler-Poincare and Lie-Poisson structure of the CH equation in ( 1.9 1 makes it clear how to generalize it to higher 



dimensions. The result is the EPDiff equation in n dimensions, given for to = by 

Sh 

dtmi = — (djirii + mjdAu 1 with velocity components - — = v? , (1-10) 

dm j 

where i,j = 1, 2, . . . , n. 

The bi-Hamiltonian property of CH. The one-dimensional CH equation ( |1.1[ ) may be written in bi-Hamiltonian 
form as 

m t = -{d x - dl) 5 ^ = - (d x {m + .) + {m + u)d x ) , (1.11) 
where the two Hamiltonians are given by 

Hi[m] = \ j mudx and Hzlm] = \ j (u 3 + uu x + 2uju 2 )dx . (1-12) 
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The integration is over the real line, for functions that decay sufficiently rapidly as |x| — > oo, and over one period, 
for periodic functions. By Magri's theorem B31 . the bi-Hamiltonian property of CH implies an infinite sequence 
of conservation laws, obtained by a recursion relation. 

The Lax pair for CH. Its bi-Hamiltonian property also implies that the CH equation < | 1 - 1 p admits a Lax pair 
representation, given by ||6l|7l 

*xx = Q + A(m + w))^, (1.13) 
% = (W-Ujttas + ^T* +7*. (1-14) 



>2A 

where u, 7 are arbitrary real constants and the eigenvalue A is independent of time. The compatibility of the Lax 



pair for constant A means that the eigenvalue equation in ( 1.13 1 is isospectral. That is, its spectrum is invariant 
under the flow of the CH equation. 

Many papers have been written to explore various features of the original single-component CH equation. A brief 
history of its exploration is recounted, for example, in |26l . See also [25] for a recent discussion of its singular 
peakon solutions. 

Plan of the paper 

Section [2] continues the introduction of our subject by discussing the extension of CH to the integrable CH(n,k) 
systems consisting of n components (momentum densities) and I < \k\ < n velocities that result from the isospec- 



tral problem in (2.1 1 and (2.2i. The CH(n,k) hierarchy may be written in a compact universal form that aids in the 
physical interpretation of its various equations as continuum flows and is reminiscent of the Virasoro structure of 
the one-component CH equation found in Il49t . Section [3] reviews the properties of the two-component CH(2,1) 



system in (3.1 )-( |3.2| ) and then discusses the CH(2,2) system, both of which are fluid systems. Section [4] provides 



two other examples of equations in the CH2 hierarchy: (i) the CH(2, - 1) system, in the position of the CH hierar- 
chy corresponding to the Dym equation in the CH hierarchy; and (ii) the CH(2,1) system with two time variables. 
Section [5] closes the paper by giving a brief summary of its main points and indicating some open problems for 
future research. 



2 A hierarchy of multi-component integrable extensions of CH 

The main feature of the Inverse Scattering Transform (1ST) for the two-component generalization of the CH equa- 
tion (CH2) is that its spectral problem is Schrodinger's equation with an 'energy dependent' potential, in which 
higher-order powers of the spectral parameter appear. For the history and development of the 1ST method with 
energy dependent potentials, one may consult BT1 1381 1531 l2l l47ll and the references therein. 

A hierarchy of multi-component generalizations of CH may be obtained by considering an extension of the Lax 



pair ( 1.13 1, ( 1.14 1 that preserves its form, but replaces its coefficients by polynomials in the scattering parameter A, 



as in [35], 

*xx = Q(x,X)^, (2-1) 
*t = -U{x,\)^x + \ 
where the potential Q(x, A) has the following energy dependence 



*t = -U(x,X)^ x + ^U x {x,X)^, (2.2) 



Q{x,X) = X n q n (x)+X n - 1 q n - l (x) + ... + Xq l (x)+ 1 A , (2.3) 
U(x,X) = Uo ( x ) + Jj-L + ...J^-L. (2.4) 
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The compatibility condition for d2T|), (|2.2|) gives the following equation, 



Qt + (d x Q + Qd x )U = -U xxx , 



(2.5) 



whose form is reminiscent of the Virasoro structure of the one-component CH equation found in [49]. 

Following this important clue to the nature of the equations in this hierarchy, we rewrite the equation geometrically 
as the Lie derivative Cjj of the 1-form density (Q dx 2 ) with respect to the vector field U as 



(d t + Cu)(Qdx 2 ) = ^(dU xx )dx, 
whose relation to continuum flows may be emphasized by rewriting it in characteristic form as 



d 



Q(x(t),t) dx(tf 



1 dx 
-{dU xx )dx{t), along — = U(x(t), t). 



(2.6) 



(2.7) 



Interpretation. The left hand side of equations (2.6 1 and (2.7 ) represents sweeping of the wave momentum density 
Q by the flow of the velocity vector field U, while the right hand side represents dispersion of the wave and forces 
that cause motion relative to the flow lines of U. 



CH(n,k) family of equations. Upon substituting the expansions in A and A -1 from (2.3 1 and (2.4), respectively, 
into equation (2.5 ), one obtains a chain of n evolution equations with k + 1 differential relations for the n + k + 1 
variables q\, qi, . . ., q n , uq, U\, . . ., Uk (n and k are arbitrary positive, or negative, integers): 



Qn—r,t 



^ {d x q n - s + q n - s d x ) 

s=max(0,r- fc) 

min(n,fc— r) 

= -(d x - dl)u r + ^2 (d x q s + q s d x )u r+s , 

s=l 

r = 0,1,..., k- 1 

o = (d x -d 3 x )u k . 



0,1,. ..,n- 1, 



(2.8) 



The differential relations in the middle equation of (2.8) contain the same (Hamiltonian) operators as in the bi- 
Hamiltonian structure < | 1 . 1 1| ). 

System (2.8) is similar to the hydrodynamic chain studied in a series of papers ll54l l47l l48ll . and to other CH 
generalizations ||2~T1 [BTl [T9l l42l [TOl . In the present paper, we shall examine the geometric structure and numer- 
ical solution behavior of several examples of equations from the CH(n,k) family of equations resulting from the 
isospectral problem (|2.1|) and (|2.2[). 



2.1 CH(n,l) system 



The simplest chain system with n coupled equations arises from the general framework (2.8) for k = 1. This 
system is denoted CH(n,l). In this case, the Lax pair (2.1 1 - (2.2 1 becomes 

Q(x,X) - ^ n ~ 2 ' xn ~ 1 
U(x,X) 



-X n p 2 + A r ' 
1 



■ n -l + ... + Xqi + 



2A 



+ u. 



(2.9) 
(2.10) 



There is one differential relation, q\ = u — u xx + const. Compatibility of the corresponding Lax pair in this case 
results in a chain of equations in the following form for p = 1,2, ... ,n, 



d t q P + {d x q p + q p d x )u - \d x q p+ i = with q n = -p 2 and q n+ i = 0. 



(2.11) 
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2.2 CH(n,2) system 
The Lax pair pTj ) - ( [22]) with 



Q = A n g n + A n ~ V-i + . . . + Agi + - , (2.12) 

U = + y + (2-13) 

generates the coupled n-component CH(n,2) system for p = 1,2, ... ,n, 

t q P + (d x q p + q p d x )u + {d x q p + q p d x )ui - \q p +2,x = 0, (2.14) 



with q n +i = q n +2 = 0. The differential relations in (2.8 1 give q\ and q2 in terms of uq and u\ in (3.5 1 and (3.6 1, 
below. The others (q3, ■ ■ ■ , q n ) are independent. 

2.3 CH(n,k) system 

The pattern continues for CH(n,k) upon including more velocities, uo,u\,... ,Uk-i, with u k = —\. The p-th 



equation in the integrable n-component system CH(n,k) is obtained from the Lax pair (2.1 1 - (2.2 1 as 

fc-l 

dtq P + ^2{d x q p+j + q p+j d x )uj - -q p+k , x = 0, (2.15) 

3=0 

with q n+1 = q n+2 = ■■■ = q n+k = 0. 
The solutions for uq,..., u k -i are updated at each time step from q\,...,qk by imposing the k+1 differential relations 



in (2.8 1 



In the remainder of the paper we will concentrate most of our attention on the two-component systems, including 
CH(2,1), referred to simply as CH2, as well as CH(2,2) and CH(2,-1). (In the CH(n,k) notation, negative values 



of k refer to positive powers of A in the expansion of U in (2. 13 1.) However, in the case n = 3 the examples 
CH(3,1) and CH(3,2) will reveal the graded structure of the Hamiltonian formulations of all of the equations in the 
CH(n,k) family. We will also show a few numerical solutions of these equations that provide crucial insight into 
their pulse-like evolutionary behavior, inviting further investigation. 

3 Examples of two-component CH systems 
3.1 Example: CH2, or CH(2,1) for n = 2, k = 1 

We denote uq = u, q\ = q and q2 = ±p 2 , and choose u\ = —1/2. In this notation, the CH(2,1) system can be 
written in the form 

qt + (d x q + qd x )u=fpp x = 0, (3.1) 
Pt + (up) x = 0, (3.2) 

where q = u — u xx + u and oj is an arbitrary constant. 

This is generally known as the CH2 system, and it has been studied extensively. The last term in the motion 
equation in CH2 has a choice of sign (=p). For the positive choice, the CH2 equation may be regarded as a model 
of shallow water waves lTT4ll29l . Its generalization to higher dimensions is immediate. 
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Spectral problem for CH2. The spectral problem for CH2 is, as in equations (2.9 1 - (2. 10 1, 

V xx = (-AV(a?) + Ag(x) + i)*. (3.3) 

This spectral problem is a type of Schrodinger equation with an 'energy dependent' potential, i.e. it is quadratic in 
terms of the spectral parameter and moreover the potential functions multiply the spectral parameter (the so-called 
weighted problem). There are some common features with Sturm-Liouville spectral problems, see for example ll38l 
|41] [53). An 'energy dependent' spectral problem also appears in the inverse scattering transform of an integrable 
generalization of the Bousinesq equation (Kaup-Bousinesq equation) fiTTl . 



Brief history of the CH2 equation. The system ( 3. 1 1-( 3.2 1 representing a two-component generalization of the 



CH equation was initially introduced in 11501 as a tri-Hamiltonian system. It was studied further by others, see, e.g., 
ll44l[8l|20l[H|29l|22l. The known applications of the CH2 model are the following. 

• In the context of shallow water theory, u can be interpreted as the horizontal fluid velocity and p is related to 
the water elevation in the first approximation tfl4"ll36ll . 

• In Vlasov plasma models, CH2 describes the closure of the kinetic moments of the single-particle probability 
distribution for geodesic motion on the symplectomorphisms lf55ll3Tll3"2"ll . 

• In the large-deformation diffeomorphic approach to image matching, the CH2 equation is summoned in a 
type of matching procedure called metamorphosis [33]. 



The same CH2 system appears as a member of the hierarchy of hydrodynamic chains studied in ll54ll . Its analytical 
properties such as well-posedness and wave breaking were studied in lfl8l l34l l56l l22l l9l and others. 
Geometrically, the original CH equation may be interpreted as governing geodesic motion for an H 1 metric that 
is invariant under the Virasoro group, as found in ll49l . Because of this property, CH is geometrically reminiscent 
of the Euler rigid body equations, which describe geodesic motion on the rotation group with respect to the metric 
supplied by the moment of inertia. The geometric interpretation of CH2 is similar: CH2 is the equation for geodesic 
motion on the semidirect-product Lie group of diffeomorphisms acting on densities, with respect to the H l on 
the horizontal velocity and the L 2 norm on the elevation. This is analogous to the finite-dimensional case of an 
ellipsoidal underwater vehicle (UWV), whose motion may be modelled as geodesies on the Euclidean group of 
rigid body rotations and translations. This finite-dimensional model of the UWV has the same Lie-Posson bracket 
as for the Hamiltonian description of the heavy top, so CH2 may also be interpreted analogously to the heavy-top 
equations. For additional discussions of geometric aspects of the CH2 system we refer to Il30ll43ll33l . 
In general, one can show that small initial data of the CH2 system develop into global solutions, while for some 
initial data wave breaking occurs |[T8l [T4l l34l l56l l22l [91 . It is interesting that only the plus sign (+) in (|3~T 



corresponds to a positively defined Hamiltonian and straightforward physical applications to shallow water waves. 



It would be interesting to know the physical interpretation of the model with the choice of the minus sign in (3.1 1, 
since this case is also integrable. 

Solutions of CH2 for dam-break initial conditions. Figure [3~Tj plots the evolution of CH2 solutions for (u, p) 



governed by equations (3.1 3.2 1 with the + sign choice in the periodic domain [—L,L] with dam-break initial 
conditions given by 

u (x, 0) = 0, p (x, 0) = 1 + tanh(x + a) — tanh(x — a) , (3.4) 

where a -C L. 

The dam-break involves a body of water of uniform depth, retained behind a barrier, in this case at x = ±a. If 
this barrier is suddenly removed at t = 0, then the water would flow downward and outward under gravity. The 
problem is to find the subsequent flow and determine the shape of the free surface. This question is addressed in 



Holm & Ivanov 



Multi-component CH equations 8 





X X 

Figure 3.1: Dam-break results for the CH(2,1) system in equations ( 3.ip.2| ) show evolution of the density p (left 



panel) and velocity u (right panel), arising from initial conditions (3.4 1 in a periodic domain. The color bars show 
positive density on the left and both positive and negative velocity on the right. The soliton solutions are seen to 
emerge symmetrically leftward and rightward after a finite time, and the evolution of both variables generates more 
and more solitons propagating in both directions as time progresses. Figures are courtesy of L. O Naraigh. 



the context of shallow-water theory, e.g., by Acheson ], and thus serves as a typical hydrodynamic problem of 



relevance for CH2 solutions with the + sign choice in (3.1 ). 

The CH2 system invites further generalizations and applications. For example, its two-time generalization is pre- 



sented in { 4.2 of this paper. 



3.2 Example: CH(2,2) forn = k = 2 

CH(2,2) designates the case of two momentum densities q\ and q 2 and two velocities uq and U\. The choice 



u 2 = —1/2 automatically solves one of the relations in (2.8 1. The other two differential relations in (2.8 1 can then 
be integrated spatially in x, to find the relationships between the momenta and the velocities, as 



qi = ui - ui jXX + u) 1: 

q 2 = «o - UQ,xx + 3u? 



u 



l,x 



2uiu 1}XX + 4wiMi + w 2 , 



(3.5) 
(3.6) 



where uj\ and lo 2 are constants of integration that depend on boundary conditions. For CH(2,2) the evolutionary 
system (|2.8[) yields equations for q\ and q 2 given by 



Qi,t + (d x qi + qid x )u + (d x q 2 + q 2 d x )u 1 = 0, 
q2,t + {d x q 2 + q 2 d x )u = 0. 



(3.7) 
(3.8) 



These equations may be solved by first updating q 2 , then q\ in (3.7 1 and (3.8 ), followed by inverting the Helmholtz 



operator twice, first for u\ in d33j) and then for uq in (3.6 1 



Dam-break equivalent problem in q 2 for CH(2,2). Figure 3.2 shows the evolution of pulses in the velocity 
variables uq (left panel) and u\ (right panel). These pulses arise from an initially localized disturbance in q 2 with a 
tanh-squared profile, for uj\ = u 2 = 0, 



q 2 (x, 0) = [1 + tanh (x + 1) — tanh (x 
m(x,0) = qi(x,0) = 0.01, 



i: 
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Figure 3.2: This figure shows the dam-break solution behavior of the CH(2,2) system (3.7 1-(3.8 ), which should be 



compared with the corresponding results in Figure 3.1 for the system of CH(2,1) equations in (3.1 H3.2 1. The dam- 
break results for the CH(2,2) system show evolution of the velocity variables no (left panel) and u\ (right panel), 
arising from initial conditions for J~q% similar to (3.4 > in a periodic domain. A "fan" of sequentially smaller, 



slower pulses is emitted rightward, after which the leading pulses overtake the slower lagging pulses in the periodic 
domain and suffer elastic collisions that show a variety of different phase shifts, apparently depending on their 
relative speeds. Each confined positive pulse in uq is accompanied by dipole (±) excitations in u\, whose polarity 
may occasionally reverse from (±, red-blue) to (=F, blue-red) and back. Figures are courtesy of J. R. Percival. 



that interacts with a constant mean flow in the velocity field u\. The constants lo\,uj2 are set to zero. The initial 
confined pulse q2(x, 0) corresponds to a confined pulse in uq that propagates steadily rightward and generates a 
structured dipole pulse in u\ that accompanies the no pulse, but may oscillate in polarity as it propagates through 
a background "fan" of smaller slower pulses. This is an interesting scenario whose dynamics will be investigated 
further elsewhere. 

3.3 A semidirect product interpretation of the example CH(2,2) 



The system ( |3.7p.8[ ) for CH(2,2) may be rewritten as 

(d t + C U0 )(q 2 dx 2 ) =0, {dt + £ U0 ){q 1 dx 2 )+C Ul (q2dx 2 ) = 1 (3.9) 
where C u denotes Lie derivative with respect to the vector field u. 



To understand system (3.9 1 better geometrically, we shall rederive it from Hamilton's principle with a Lagrangian 
defined on the semidirect product Lie algebra of vector fields 1(uq, u\) : jCo(K)(§)Xi(K) — > K. Here no £ Xo(M) 
and u\ G X\ (K) are right-invariant smooth vector fields on the real line R. Vector fields in Xo (K) act on themselves 
and on Xi(M) by vector cross product, while the action of Xi(IR) on itself is assumed to be by simple vector 
addition. That is, vector fields in are advected quantities, see, e.g., E81 . 

Definitions: Semi-direct product. The semidirect product Lie algebra action Xo(R)(§)Xi(K) is defined by 

[(x ,Xi),(y ,n)] = ([x ,y ],[Xo,yi] + [^ y o]), (3.io) 

where [ • , • ] is the commutator of vector fields in natural notation. This commutator defines the adjoint action 

ad (X0)Xl) (y ,yi) = -l(x ,x 1 ), (y ,yi)] = (ad Xo y , ad^y + ad Xl y ) Q.ii) 
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and the L 2 pairing with (a, f3) G Xq(S)X* yields the coadjoint action 



ad (x ,x 1 )(«'/ 3 )' ( Y o, Y i)) = ((">/5)> ad (XoiXl) (y ,y"i 
: ((a, /3), (adxo^b, adx ^i + adj^Yb)) 
■- (a, &d Xo Yo) + {(3, adxo^i + adj^^o) 
: (ad^a + ad^/3, Y ) + (ad^/3, Yi) 
^(ad^a + ad^Aad^), (Y ,Y{)) . 



(3.12) 



See J5J] for more background and derivations of the formulas for the adjoint and coadjoint actions of the semidirect 
product Lie group Diffo(M)(s)Diffi(M) and its Lie algebra of right-invariant vector fields Xo(|)Xi. 

We shall rederive system ( |3.9[ ) from Hamilton's principle by using the Euler-Poincare theory, as reviewed for 
continuum mechanics, e.g., in [28]. We then pass to its Hamiltonian formulation in terms of a Lie-Poisson bracket 
by performing a Legendre transformation. In the Euler-Poincare framework, we have the following. 

Theorem 1 (Euler-Poincare formulation of CH(2,2)). 

Hamilton's principle 5S = with S = f 1(uq,u\) dt yields the CH(2,2) system in equations (3.9) for the La- 
grangian 



l(u ,ui) 



dx, 



UoUl + Uo,z«l,z + ^>l u + LU 2 Ui + 2ui (u-l + u lr 
for constrained variations ofuo and u\ of the semidirect product form in \3.10 ), consisting of 

5u = d t £ + [uq, Co] = d t £o - ad M0 £ , 
5m = d t (,i + [«o,£i] + [ui,£o] = dt£i - ad Uo £i - ad Ul £ - 

Proof. By direct calculation, Hamilton's principle with the Lagrangian ( 3.13| > implies 

= 5S = 



dt 



(3.13) 



(3.14) 
(3.15) 



(qi, 5u ) + (q 2 , Sui) dt 
J (qi, dtto ~ ad Mo ^ ) + (<?2, dtii - ad Uf) £i - ad Ul £ 
- J {dm + ad « <7i + ad* x g 2 , £0} + {d t q 2 + ad£ g 2 , £1) dt , 



where q\ and q 2 are given in ( |3.5| ) and ( |3.6| >, respectively, and we have used the semidirect product Lie algebra 
action defined in (3.10). The coadjoint operation ad* is defined using the I? pairing as in, e.g., 



(<?2, ad M0 £i) = (ad* g 2 , 6/ 
which is a repeated pattern in the system (|3.9[). 



{£uoQ2, £1) = ((9x42 + q2d x )u , £1 



(3.16) 



□ 



Legendre transform to the CH(2,2) Hamiltonian formulation. The Legendre transform for CH(2,2) is given 
by 

<5/i 8h 51 

h(qi,q2) = ({qi,q2), ( u o,ui)) - l{u 1 ,u 2 ) , Jf n =u ' Su = 5^' ^ ' 

where ( • , • ) denotes 1? pairing on the real line. The corresponding variations yield 

5h = ((u ,ui), (Sqt, Sq 2 )) + /(qi~ q2 - , ($u , 5ut)\ . (3.18) 

Hence the pairs (qi,q 2 ) and (uq,u\) are dual variables with respect to the Legendre transform of the CH(2,2) 
Lagrangian in ( |3.13[ ). Thus, the differential relations in the Lax pair formulation of CH(2,2) appearing in (3.7 1 and 
(|3.8[) Legendre transform into dual momenta for the Hamiltonian formulation. 
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Lie-Poisson Hamiltonian form of CH(2,2). The CH(2,2) system in ( |3.7p.8[ ) or ( |3.9[ ) may now be cast into 
Lie-Poisson Hamiltonian form, as 



d t qi 









d x qi + qid x d x q 2 + q%d x 
d x q2 + q2d x 



5h/5q\ 
Sh/ Sq2 



u 



(3.19) 



The Hamiltonian operator yields the Lie-Poisson bracket defined on the dual to the semidirect product Lie algebra 



of vector fields £o(M)(S)£i(M). The Lie algebra action for the semidirect product is defined by (3.10l and dual 
coordinates are qi 6 Xq(R) and q 2 £ 3t|(M). 



Remark 1. The Lie-Poisson Hamiltonian form (3.19) 



d t (q2dx 2 



-C uo (q 2 dx 1 



on its 



may be interpreted as saying that the 1-form density (q2dx 2 ) evolves in time t by the action of Diffo( 
initial conditions. That is, 

— (q 2 dx 2 ) = along — = u (t,x(t)) . 

Quantities that evolve this way to remain invariant along the characteristic paths of a flow velocity in ideal fluid 
mechanics are said to be advected, or frozen into the flow. Equation ( |3.9[ > shows that the process is not passive, 
though; because the dynamics of q\ is affected by a force depending on q 2 and the corresponding velocities that are 
obtained from the differential relations in \2>.l) and ( |3.8| ). 

The other Hamiltonian structure for this example would be interesting to know. However, knowing it is not nec- 
essary for the sake of generating its integrable hierarchy, because we already have its isospectral problem and Lax 
pair. 



Analogy with rotating tops. If the CH(2,2) problem specified here for a right-invariant Lagrangian on X®X had 
been expressed instead on so(3)(S)so(3) for a left-invariant Lagrangian, the result would have been interpreted as 
the dynamics of a rotating top in a potential force field, as discussed in [3 ]. The dynamics in this case is expressible 
in Hamiltonian form as 



qi 




qi x q 2 x 




5h/5qx = u o 






q 2 x 




5h/Sq2 = ui 



(3.20) 



for angular momenta (qi,q2) G R 3 x R 3 and their corresponding angular velocities (uo,ui) G R 3 x R 3 and 
Hamiltonian h = ^(qi • uo + q 2 • ui). This Hamiltonian matrix defines a Lie-Poisson bracket on the dual of the 
semidirect-product Lie algebra R 3 (s)Rq in which the first R 3 acts on itself and on the second, Rjj, by vector cross 
product, while the action of the second Rjj on itself is by simple vector addition. 



4 Other examples of two-component CH generalizations 

In this section we discuss two other examples of equations in the integrable CH2 hierarchy. These are: (i) the 
equation denoted CH(2,-1) in the CH2 hierachy that occupies the "Dym position" in the KdV hierachy; and (ii) the 
CH(2,1) equations with two time variables. 



4.1 The CH2 Dym equation or CH(2, -1) 



The CH(2,2) system ( |3.5| ) - ( |3.8| ) sits in a hierarchy of integrable equations that share the same spectral problem 
(3.3), as well as the main representative CH(2,1). There are other members of this hierarchy, for which U(x, A) 
contains positive powers of A. A simple interesting example of this kind has 



Q = - X 2 p 2 + Xq + l/4 and U = \u , 



(4.1) 
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Extending our notation, we can denote the resulting system as CH(2,-1){^] Substituting (4.1) into (2.5 1 leads to 
equations denoted as CH(2, -1): 



Qt 



Pt + 



;i-^ 2 ) 



0. 



(4.2) 
(4.3) 



where we have used up = K, obtained from the differential relation arising in the A 3 term, and have set the constant 
value K 



-2. Some of its solution behavior is shown in Figure 4.1 





Figure 4.1: Results are shown for the evolution of the CH(2,-1) system (4.2 4.3) with e = — 1 for p (left) and q 
(right), arising in a periodic domain of length L = 80 from initial conditions that represent a dam-break q(x, 0) = 
tanh((x — L\)/a) — tanh((x — L2)/a), p(x, 0) = 1 with a = 1, L\ = L/3, L2 = 2L/3. Soliton solutions are 
seen to emerge and propagate in both directions. The head-on collision process produces a slight refraction of the 
soliton trajectories, unlike the CH(2,1) case. Figures are courtesy of V. Putkaradze. 

Remarks. 

Here are a few remarks about the CH(2,-1) system in equations (|4.2|) and (|4.3[): 



• This coupled nonlinear system is at the position in the CH2 hierarchy that corresponds to the modified Dym 
equation, first introduced as a tri-Hamiltonian system in ll50l . 



The CH(2,-1) equations combine to produce the nonlinear wave equation 



(4.4) 



Linearizing this equation around q = and p = 1 yields the dispersion relation for a plane wave ex.p(i(kx 
cot) with wave number k and frequency co as 



oo 2 (k) 



[l + k 2 )k 2 . 



Accordingly, the phase speed of the linearized plane waves is oj/k = y/l + k 2 , so the higher wave numbers 
travel faster. This type of dispersion relation is not unfamiliar: it is the same as for time-dependent Euler- 
Bernoulli theory for an elastic beam with with both bending and vibration response [52]. 



3 It is possible to extend the expansions in powers of A in 1 2.3 1, i2.4i in both directions for both Q and U. The equation chosen for 
analysis here is only a single step in this extension. 
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The CH(2,-1) travelling wave solutions v(£) = l/p(£) with £ = x — ct for c > conserve the energy E 
given by 



2E={v'f-(l-l) 2 -tv + J?, 
with integration constants / and J defined by 

C n it" 

7 = gt> and J = v — v + cq . 

v 

For the travelling wave £7 = 0, and when / = = J, as well, then the solution is given by 

cp(0 = \/sech(2(£ - &)) , Co = constant and qr(f) = c/5 3 (£). 
This is a confined travelling wave pulse in both p and q. 

The coupled CH(2,-1) system ( 4.2||4.3 1 may also be written in Hamiltonian form, as 







~d x 




—q/p 2 = 5h/6p 


e 1. 




. o d x - dl_ 




l/p = 5h/8q 



d t 



This Hamiltonian operator yields the Poisson bracket dynamics, 



with h 



(q/p) dx . 



d£ 
~dt 



('■->-/(^^(*-*t)*- 



The energy conservation law may be expressed in conservative form as 

2 



0, ( J ] + <),,. 



1 (± 

2 {pi 



+ 



1 

2? 



p P . 



0. 



The CH(2,-1) equations (4.2) and (4.3 1 may also be written in n dimensions as 

« + v-(£) = o, 

q,-v((l-V 2 )I) = 0, 

for p < R and q G M n . 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(4.11) 
(4.12) 



The consistency among the CH(n,k) equations can be demonstrated b y com bining them. The choice, Q = 
e\\ 2 p 2 + \q + \ and U = — ht + u + 2e2\/p with ei^ = ±1 in (2.5 1, for example, produces another 
nonlinear integrable system]^] 



ei (jh + (up)x) + £2 



% + + qd x )u - eipp x + e 2 (l - d\ 



P 2 , 



P 



o, 

0. 



(4.13) 
(4.14) 



with differential relation q x = u x — u xxx . This system reduces to CH2 for e 2 = 0. 



One may compare this with Q = - A V + Aq + 1/4 and U = Xu for CH(2,-1) in J41 
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4.2 Equations in the CH2 hierarchy with two time variables 



There is also an integrable CH2 system with two 'time' variables (t and y). In particular, consider the systerrj^The 
two-time CH2 system with m = U x — U xxx is 



m t + 2U yx m + (U y + j)m x + pp y 

This system can be written equivalently in a hydrodynamic form as 

{m/p 2 ) t + (U y + j){m/p 2 ) x = - 
Pt+({U y + l)p 



0, 
0. 



P Py, 



0. 



(4.15) 
(4.16) 



(4.17) 
(4.18) 



which shows that it has only one characteristic velocity, dx/dt = (U y + 7). 

The two-time CH2 system can also be written as the compatibility condition for the following linear system (Lax 
pair) with a constant spectral parameter (: 



2C 



(Uy+j)^ X + -U yX ^. 



(4.19) 
(4.20) 



The first equation in this system is the spectral problem (3.3 1 of the CH2 hierarchy. The second equation introduces 
the other 'time' derivative, with respect to y. The system (4.15 i,( 4.16| > appears on setting 



then using (4. 19 1,( 4.20 > to eliminate higher derivatives and assuming Q = = Q. 
Remarks. 



Perhaps not unexpectedly, the corresponding modification of the linear system (2.1 1, ( |2.2| ) yields a two-time 
version of the entire CH(n,k) hierarchy in (|2.8|>. 



• The integrable system of two-time CH2 equations (4. 15 1,(4. 16 1 reduces to CH2 for x = y and u = U x . 

• Likewise, the special case 7 = = to with initial condition p = admits iV-peakon solutions, 



A' 



m{x, t, y) = ^2p a {t, y) 8{x - q a (t, y)) 



(4.21) 



0=1 



with two 'time' variables (t and y). 



The amplitudes of m and p as functions of x and y in Figure |4.3| show modulations along the crest of the 
two-dimensional solitons as they form under the two-time CH2 dynamics. . 



5 A two-time version of the CH equation has been considered previously in 1371 . 
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0.025 




1.3 






2.6 





Figure 4.2: An initially sinusoidal wave train (plotted with m above and p below at five values of y shown in 
colou r s) con centrates into steeper, larger nonlinear waves under the dynamics of the two-time CH2 equations 
(4.15 >,( 4. 16 >. Figures are courtesy of L. O Naraigh. 



Graded Lie algebra structure for the CH(n,k) chain with n > 3 and k > 0. A new feature of the CH(n,k) 
chain may be recognized for n > 3 and k > 0. In that case, the Hamiltonian operator reveals its character as the 
Lie-Poisson operator defined on the dual space of a graded Lie algebra. For example, in the case of CH(3,1) the 



system in the first line of (2.8 1 may be written in Lie-Poisson Hamiltonian form, as 



~d t q{ 




d t q2 




Pt<&_ 





d x q\ + qid x d x q 2 + qzdx d x q 3 + q 3 d x 
d x q2 + qA d x q 3 + q 3 d x 
d x q 3 + qzdx 



8h/5qi - 
Sh/5q 2 = 
5h/5q 3 



u 
_i 

2 

= 



(4.22) 



For q 3 = —p 2 this generalizes the CH(2,1) equation to three components. In hindsight, we see that we could have 
written the CH(2,1) equation in the same Lie-Poisson Hamiltonian form, as 



dtqi 




d x qi + qid x 


d x q2 + qzd x 




5h/5q\ = uq 


d t q2_ 




d x q2 + qA 







5h/5q 2 = — | 



(4.23) 
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0.8 

0.6 





Figure 4.3: An initially sinusoidal wave train in the (x, y) plane (whose amplitudes are shown by colour bars) 
with m above and p below, steepens and modulates along th e cres t s as it grows into a sequence of nonlinear wave 
packets under the dynamics of the two-time CH2 equations (4. 15 1,(4. 16 1. Figures are courtesy of L. O Naraigh. 



for the Hamiltonian h(qi, q 2 ) = | / <7i(l — d%) 1 qi — q2 dx. Likewise, in the case of CH(3,2) the system becomes, 



~d t qi 




d t q2 




p t q 3 _ 





x qi + qid x d x q 2 + q 2 d x d x q 3 + q 3 d x 
d x q2 + q2d x d x q 3 + q 3 d x 
d x q 3 + q 3 d x 



Sh/Sqi 
Sh/5q 2 
5h/5q 3 - 



Ul 



(4.24) 



qi,q2,q3 



This grading of the Hamiltonian operator according to the weight n reveals the new feature. For CH(3,k) with 
k > 0, the Hamiltonian operator in ( |4.24 i defines a Poisson bracket on the dual of the Lie algebra of three- 
component vector fields (Xi, X 2 , X 3 ) G £i x £2 x £3 defined by their graded commutation relation 

[(X 1 ,X 2 ,X 3 ),(Y 1 ,Y 2 ,Y 3 )} = ([X^YxUXx,^] + [X 2 ,Fi], [X h Y 3 ] + [X 2 ,Y 2 ] + [X 3 ,Yi]). (4.25) 

Dual coordinates are q\ G 3£|(]R), q 2 G £|(K) and q 3 G £|(R). The Lie-Poisson bracket above may now be 
written as 

'SF 5F 5F\ f5H 5H 5H" ' 
^Sqi' 5q 2 5q 3 J ' \oq : <u !2 oq 3 _ 

in terms of the graded commutation relation (4.25 1 and the L 2 pairing ( • 
its dual. 

Remarks. 

• A Lie-Poisson bracket defined on the dual of the same graded Lie algebra also appears in plasma theory, at third 
order in the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy of equations [46]. 

• From the graded Lie-algebra action ( |4.25 ) and the Lax pair in ( |2.1| ) - ( |2.2| ), it is clear how to extend the pattern to 
higher order and thus include the more deeply nested systems in the CH(n,k) chain. 



(4.26) 

between the graded Lie algebra and 



• The semidirect product action in (3.10) is also an instance of the graded Lie-algebra action. 

• This type of weighted Lie-Poisson bracket is also encountered in the Hamiltonian formulation of dynamics in the 
BBGKY hierarchy for ideal plasma physics [46]. 

• A graded extension of ( 3.20[ ) analogous to that for the Hamiltonian operator in (4.24 1 is also available for coupling 
with additional so(3)* angular momenta. 



• The properties of the full Hamiltonian structure for the CH(n,k) chain in ( 2.8 1 consisting of n evolution equations 
with \k\ differential relations that allow k < will be discussed elsewhere. (For k < 0, the grading runs in the 
'opposite direction' in a certain sense.) 



5 Conclusion 



Main results of the paper. We first formulated the Lax pair consisting of the energy-dependent isospectral prob- 
lem and evolution equation (2.1 1, (2.2 1, whose compatibility yields the integrable family of CH(n,k) systems ( 2.15[ ) 
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with n components (momenta) and 1 < \k\ < n velocities. After looking at several examples among the CH(n,k) 
multi-component equations, we investigated some of the other equations of the CH2 hierarchy and found geomet- 
rical comparisons with systems of coupled spinning tops, as well as fluids because of the semidirect-product nature 
of their Lie-Poisson brackets. In particular, the integrable CH(2,2) equations in (3.9 1 with Hamiltonian matrix in 
equation ( 3. 19 1 were seen to be analogous to the finite-dimensional equations for a spinning top in a potential force 
field J3j|23l|24l, whose Lie-Posson bracket is dual to the semidirect-product Lie algebra so(3)(|)so(3)o, in which 
the second entry is treated simply as a vector space, as discussed in j ]3.3| 

Section ^provided additional examples of other integrable equations in the CH2 hierarchy, such as the CH(2, -1) 
system and the two-time version of the CH(2,1) equations. The CH(n,k) systems arising from ( |2.1| )- p3] > with 
negative values of k were found to show quite different Hamiltonian structures from their corresponding systems 
with positive values of k. 



Properties of the CH(n,k) systems. Several properties of the CH(n,k) systems were identified in the course of 
this work. These included the following. 



(i) The differential relations in the middle equation of (2.8 1 involved both of the compatible Poisson operators 
in the biHamiltonian structure ( |1.11| >. In both CH(2,1) and CH(2,2) the differential relations defined the 
momenta dual to the velocity vector fields. 



(ii) The Hamiltonian structure for the CH(n,k) chain with k > in (2.8 1 consisting of n evolution equations and 
k + 1 differential relations possesses a graded Lie algebra structure, which became evident for n > 3. In 



hindsight, looking at (3.19), the semidirect-product Lie algebra structure for n = 2 in (4.23 1 could have 



already been understood as being graded. For k < again both of the compatible Poison operators appear. 

(hi) The sample numerical simulations shown for the CH(2, 1), CH(2,2) and CH(2,-1) equations revealed challeng- 
ing properties for future investigation, such as different types of collision behavior. In the case of CH(2,1) 
with two times, the simulations also revealed modulation of the waves along their crests during the formation 
of the soliton trains. These sample numerical simulations provided insight into the fascinating pulse-like 
solution behavior of the CH(n,k) equations and invited further investigation. 



Future challenges. The CH(n,k) family of integrable partial differential equations (PDE) discussed here offer 
many interesting challenges for future research, particularly in determining and analyzing their solution behav- 
ior and possible physical applications. Besides CH(2,1) which may be interpreted as a shallow water system the 
CH(2,2) equations in ^J3} the CH(2,-1) system in ^4.1| and the two-time CH2 equations in {4.2 all offer new chal- 
lenges for physical interpretation and mathematical analysis. For example, one may expect the continuing interest 
in wave-breaking analysis for CH and CH2 to extend also to the other integrable PDE in the rest of the CH(n,k) 
family, including, e.g., the CH(2,-2) system, which was not discussed here. The numerical simulation of these in- 
tegrable PDE, and the formulation and analysis of their discrete versions can also be expected to attract attention in 
future endeavors. Finally, the multi-dimensional extensions and deeper geometrical aspects of these new integrable 
PDE also pose interesting challenges for future research. 
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